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"Ti • Abstract 

^ • 

f-\ . Studies of three-dimensional and two-dimensional condensed phases have shown that many- 

O ' body interactions contribute ~ 10% to the equations of state of noble gases. This paper assesses the 

importance of three-body triple dipole interactions for quasi-one-dimensional phases of He, Ne, H2 , 

Ar, Kr and Xe confined within interstitial channels or on the external surfaces of nanotube bundles. 

We find the substrate-mediated contribution to be substantial: for interstitial H2 the well depth of 

the effective pair potential is reduced to approximately one half of its value in free space. 

*yv ' We carry out ab initio calculations on linear and equilateral triangular configurations of (H2)3 and 

f^ I find that overlap interactions do not greatly change the DDD interaction in the linear configuration 

J~^ , when the lattice spacing is greater than about 3 A. However, the DDD interaction alone is clearly 

C^ ' insufficient for the triangular configurations studied. 
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1 Introduction 



Adsorbates within carbon nanotubes and bundles of nanotubes have been the focus of much recent 
attention (see for example refs. B-Pl). In particular, novel quasi-one dimensional adsorbate phases 
have been predicted. The closest to true one-dimensional behavior is likely to be realized by small atoms 
Q \ or molecules (He, Ne, H2) in the "interstitial channel" (IC) phase, and the "groove" phase. For the IC 

O ' phase, the molecules adsorb in channels bounded by three contiguous nanotubes within bundles. The 

groove phase comprises adsorbate in the grooves between adjacent tubes on the perimeter of a nanotube 
bundle. In both of these phases the molecules are confined transversely but can move more or less freely 
along the lengths of the nanotubes. 
jrt ' To date, in exploring adsorption on and in nanotubes, relatively little attention has been paid to the 

detailed nature of the interactions among adsorbate molecules. Indeed, most studies have assumed that 
the pairwise interaction does not differ significantly from its free-space form. In this paper we show that, 
for the IC and groove systems considered, the proximity of a polarizable medium, the nanotube walls, will 
profoundly alter the effective pair potential between two adsorbate molecules, to the extent of producing 
qualitative changes in expected behavior. 

In two- and three-dimensional systems many-body effects have been extensively studied P, fiO, |ll| . 
Barker |TQ], in particular, showed that equilibrium properties of condensed rare gas systems could be 
obtained accurately using potentials in which the pair potentials were obtained from gas-phase data, and 
the higher-order contributions were purely of Axilrod-Teller-Muto or "triple-dipole" (DDD) type, even 
though this is exact only asymptotically for large molecular separations. Barker argued that extensive 
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cancellation must take place among other many-body interactions, at least one of which is known to be 
comparable in magnitude to the DDD interaction. 

Analogously in two dimensions, adsorbed monolayer films manifest two types of many-body con- 
tributions that are relevant to quasi-one-dimensional nanotube phases. One is that involving three 
adatoms. The other is a modification of the adatom pair potential by the substrate acting as a third 
body. McLachlan derived an asymptotic expression, valid at large adatom separation, for this substrate 
mediation energy [n2J. As in the three-dimensional case, experimental data at condensed phase densities 
in adsorbate films are consistent with the asymptotic expression [|l^ [Ij]. However, the experimental data 
for two-dimensional films are too limited to provide a thorough test of the adequacy of McLachlan's in- 
teraction. We note that a third three-body energy in films involves two substrate atoms and an adatom. 
This oft-neglected interaction can make a substantial correction to the gas-surface adsorption potential. 

While extensive cancellation occurs as discussed above in three dimensions, among many-body inter- 
actions other than the DDD three-body interaction, this cannot be assumed a priori in one dimension 
because of the strong angular dependence typical of many-body interactions. We have therefore extended 
a recent accurate ab initio calculation of the potential energy surface of H2-H2 [ fi5[ to linear equispaced 
(132)3 and to (132)3 arranged in equilateral triangles. While significant deviations from the DDD inter- 
action are found for the equilateral triangular configurations, the deviations from DDD for linear (£[2)3 
are important only for separations small compared to the equilibrium separation of H2~H2. 

In this paper we consider a linear array of rare gas atoms or of H2, henceforth denoted A, each 
separated from its nearest neighbor by a lattice constant a. In Section 2 we evaluate the net DDD 
contribution arising from interactions among all sets of triplets AAA. The resulting three-body energy 
per atom, which we call E)^j^^, turns out to be negative and of order 1% of the total pair interaction 
energy. We next evaluate an energy per pair V^^q, due to a summation of the DDD interactions between 
two A particles and all C atoms. This energy is found to be positive and of large magnitude, equal 
approximately to half of the interaction between the atoms in free space. When Vj^J^q is added to 
the free space interaction, there results an effective pair potential which has a well depth significantly 
smaller than the free space well depth. This reduced attraction implies that collective phenomena|^, 
such as condensation, ought to occur at much lower temperatures than would be predicted if only the 
free-space interactions were taken into account, and substrate mediation were ignored. Qualitatively 
similar behavior is derived as a consequense of the DDD interaction for He, Ne, H2, Ar, Kr and Xe 
atoms/molecules adsorbed into the grooves of a nanotube bundle. In Section 3 we discuss the ah initio 
calculation carried out on (H2 )3 in linear and triangular arrays. 

While this paper evaluates electronic substrate mediation of adsorbate interactions, a substrate phonon 
mediation mechanism also exists. Such phonon mediation effects on adatom interactions on planar 
surfaces are reviewed by Gottlieb and Bruch |lj]. While no such evaluation has been performed for 
atoms or molecules in nanotubes, the analogous effect of nanotube phonons on the electron-electron 
interaction has recently been assessed by Woods and Mahan p6| . 

2 Adsorbates in IC and Groove Phases 

We begin by considering particles A confined to the axis of an IC phase and separated by a lattice constant 
a. There are two relevant three-body energies which we consider. The first is the net DDD energy arising 
from interactions among all sets of triplets on the line. The second is the total DDD interaction involving 
an adatom pair and all possible C atoms. The potential energy $ corresponding to a linear array of N 
particles specified by their position vectors rl, ...,rjv is given by 

<l>(rl, . . . , rKr) - ^ (7(f](rl, r,) + ^ uf^M^r.^r^) + ^ ufj^^kr,, T,,r-i) , (1) 

i<3 i<i<k i<j<k 

where C/*^^^ and C/*^'^^ are pair and DDD triplet potential functions and four-body (and higher order) 
interactions are ignored. For purposes of estimating the total pair energy we will temporarily assume 

^ In this paper we do not consider the implications of substrate-mediated adsorbate interactions for condensation involving 
interactions among different IC or groove phases in a nanotube bundle (see e.g. Ref.|p|). 



that the pair potential C/)j^^ (r^ , nj ) is isotropic and of Lennard- Jones 12-6 form with distance parameter 
a and energy parameter e. The total pair interaction energy per particle A can be written as 



Ef. = 4^ = 46 



^AA 



N 



©"«-)-©'««) 



(2) 



in the limit N ^ oo, where C is the Riemann zeta function. The right-hand side of Eq. (0) has a 
minimum at a/a ~ 1.122. Evaluating Eq. (H) at that minimum we obtain: 



.(2) 

E)^/^ = ^ = -1.035e = -0.509^ , (3) 



^i^)YAA=-lM5e= -0.509^ 

^^ N a^ 



where Cg — Aecr^. Next we assess the contribution of the DDD interactions U^j^^j^ to V^^- A simple 
result for U\J^y^ holds for isotropic oscillators, at large separations: 

(3) 3cos(6>,)cos(6>j-)cos(6>fc) + l 

UaaA = ^AAA 3 ^3 ^3 ' (4) 

ij ik jk 

where vaaa is the triple dipole dispersion energy coefficient, rij,rik,rjk are the interatomic distances 
in a given triplet, while 0i,6j and 6^. are the internal angles of the triangle formed by the atoms i,j 
and k. The sign of the triple dipole energy depends on the geometry of the triangle formed by the 
three atoms. Indeed, the DDD term is positive for acute triangles, while for most obtuse triangles it 

(3) 

is negative. In particular, C/4^^ is negative for the linear configuration of adsorbate particles, since 
cos{9i) cos{9j) cos{9f^) ~ —1 in that case. Then, the net DDD potential becomes: 

Vaaa^- Y] ^!^ddi_. (5) 

i<j<k (^vnkTjk)^ 

Summing over all possible triplets we obtain the net DDD energy per atom: 



j-,(3) _ _2l^AAA 
^AAA — „9 
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^C(9) + 22^ 2^ ^3(^ + ^_l)3(„ + 2m-l)3 
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(6) 



^(3)^_a2^_ (7) 

The coefficient vaaa can be estimated from a Drude model of isotropic oscillators from the coefficient 
Cg : Vaaa = \C^olai where oca is the static polarizability of the adatom. Then the ratio of the DDD 
energy to the pair interaction energy is: 

^AAA - ^ . a 



= 0.4-. (8) 

This ratio is very small (0.8%, 0.3% and 0.5% for II2, He and Ne, respectively, using values of a and a 
taken from |11[|). 

We next evaluate and assess the net DDD interaction involving the adatom pair and all possible C 
atoms. The effective pair potential has the form: 

Kff^(^i'^-) = vfAin^i) + yflc^nyo) ■ (9) 

('2^1 {1\ (1\ 

For a linear configuration of adsorbate particles: C/j;^^ (f^ , rj ) = U\j^(\zi — Zj\) and V^^g (?^, ''j) = 
= V^jj (|zi — Zjl), where z axis is along the axis of the IC. The DDD energy Vaaq is given by: 

yAAc{^^yi) = E VfAC^^^y.y^) ■ (10) 

k 



Here f% is the position of the fc*'' carbon atom along the nanotube surface and f/y^^c" ^^ ^^^ DDD interaction 
with the adatoms at fl and fj. In the approximation of isotropic oscillators, the expression for the DDD 
energy is analogous to Eq.(0), except that the strength coefficient is: 

ia\acEcEA{Ec + 2EA) ,„, 

"■^^^ ^{Ea + EcY ■ ^''^ 

Here the energies Ea and Ec are characteristic energies of the adatom and C atom, respectively, and ac 



is the static polarizability of a C atom M, 18 



For simplicity, we perform the summation of gas-surface DDD interactions (10) by smearing out the 
carbon atoms in the substrate. This continuous approximation introduces an error, but the qualitative 
trends are expected to be accurate in general. Under this assumption, the net DDD contribution from 
three nanotubes of radius R takes the form: 

(3) _ \2-dvAAcM{x- v) 

where i? = 0.38 A~^ is the surface density of C atoms and z is the distance between adatoms. Here, 
X = z/R and y = d/R, where d is the distance from the axis of the cylinder to the axis of the IC. The 
integral M{x] y) is a dimensionless integral over a cylindrical surface. Results for V^^<^ as a function of 
adatom lateral separation, in the case of quasi-ld phases of He, Ne and H2 adsorbates, were calculated 
from (|lj) for nanotubes of radius 6.9 A and d = 9.815 A. 

The pair potentials employed here are those of Silvera and Goldman ||l^ (H2); Aziz, Nain et al. [g^ 
(He) and of Aziz, Meath et al. |£1| (Ne). The values of DDD strength coefhcient ( [ll| ) were calculated 

from Ref.O, 18|. Results of our calculations are given in Fig. 1, Fig. 2 and Fig. 3 and in Table 1. In 

(2) 
Table 1, e and Zm are the well depth and the location of the potential minimum for V^g ; eo and zom are 

the corresponding quantities for the free space potential; S is the percent decrease in the well depth due 

(3) (2) 

to Vj^^fj and p is the sign-reversal point of V^j^ . 

In all three cases we observe a large repulsive DDD effect: the reduction in the well depth of the pair 
interaction is 54%, 28% and 25% for H2, He and Ne, respectively (Table 1). An interesting feature of the 

effective potential is the presence of a critical value of adatom separation for H2 and He adsorbates, beyond 

(2) ■ • . ' % 

which V^g becomes repulsive for a certain range of separation values. These points are p ~ 6.3 A and 

7.4 A for H2 and He, respectively (Table 1), with corresponding ranges of repulsive V^g : Az=(6.3— 15.1 
A) and A2=(7.4 — 11.7 A). On the other hand, an investigation of the asymptotic behavior of M(x; y) 

reveals that the DDD energy becomes negative for very large interparticle separations ( z > 9R ). We 

(2) 
note that a reversal in sign for V^g is not found for Ne. To elucidate this fact we investigate the relative 



strength of DDD interaction at zom- The expression ( |12D can be used to yield: 

YMd^^kM{xo^;y)['^^^£^], (13) 

Co L'e 

where eo — Ce/2z^ is the free space well depth and k is independent of the adsorbate species. The integral 
M{xom', y) has only a weak dependence on zom in the region of interest (3 — 3.5 A). Then the relatively 
small DDD effect for Ne is explained by the bracketed term involving vaac , Ce and zom in (^j) • Its value 
is nearly twice as large for H2 as for Ne ( 34.1 A^ vs. 17.7 A^ ). 

The present results imply that the substrate mediation effect is approximately three times as large 
for the IC geometry than for adsorbed monolayer films pTl. Two factors are responsible. One is that the 
equilibrium distance to the substrate (2.9 A) is ~ 5% smaller than on a planar surface, yielding a larger 
substrate contribution. The more important factor is that three adjacent surfaces contribute in this IC 
geometry. 

One significant consequence of the reduced attraction between the particles within an IC is a qualita- 
tive change in the predicted equation of state for this Id system. It has been shown |g^ that the ground 
state of a Id bose system is a liquid if and only if a bound dimer exists. For He, the dimer is bound by 
just a few milliKelvin, if the free space potential is used gj, ^ |2^. There is no doubt therefore that 
this binding disappears when the weakened attraction is substituted. For H2, the dimer binding energy 



computed with the Silvera-Goldman potential is large (0.384 meV) and the ground state of the system is 
a liquid with cohesive energy ~ 0.4 meV per molecule |g^. When the new potential appropriate to the 
IC is used, we find no bound states. ( For D2, the dimer is very weakly bound, by 0.0513 meV ). If our 
calculations are correct, they imply that the ground states of interstitial "^He and para-H2 are quantum 
gases. Such a system is of considerable fundamental interest. 

In the nanotube bundle geometry, adsorption can also occur on the outer surface of a bundle. The 
most favorable site energetically is the groove formed at the intersection of two external tubes, which 
has been studied by several groups both experimentally and theoretically. In contrast to the IC case, 
larger particles ( such as Ar, Kr and Xe ) can fit nicely into the grooves M. We have computed the 
three-body effect of the substrate, using the DDD expression as in the IC's. The resulting effect on the 
pair potential for H2, He, Ne, Ar, Kr and Xe is reported in Table 2. The substrate effect of reducing the 
well depth is smaller than in the IC case. The reasons are twofold. One is simply that only two nanotubes 
contribute here, instead of three in the IC case. The second reason is that the distance to the tubes is 
somewhat larger than in the IC case; in these grooves, the distance is very close to the a parameter in 
the gas-carbon pair interaction. (It would be exactly a if the tube radius were infinite). These differences 
parallel those used above to distinguish the IC case from the graphite case. For this reason, the groove 
results are intermediate between those from the latter two geometries. 

3 Three-body interactions for the H2 trimer in free space 

We have calculated the nonadditive, three-body interactions for the linear and equilateral triangle ar- 
rangements of three hydrogen molecules from first principles. The H2 bond length is kept fixed at 0.7668 
A, consistent with the rigid rotor approximation. With the molecular centers-of-mass forming the vertices 
of a triangle, the intramolecular orientation of each molecule is described relative to three axes centered 
at each vertex. We have explored the effects of molecular orientation on the three-body interactions by 
orienting each molecule in one of three ways. Each molecule is oriented with its molecular axis perpendic- 
ular to the plane of the triangle, bisecting a corner of the triangle, or perpendicular to the line bisecting 
a corner. In total, there are 27 possible combinations of trimer "orientations" . Due to symmetry, it is 
necessary to calculate only 10 non-degenerate orientations for each geometry. The closest center-of-mass 
to center-of-mass distance was stepped by 0.2 angstrom from 2.0 to 6.0 angstrom for the linear geometry 
and from 3.0 to 6.0 angstrom for the equilateral triangle arrangement. 

Convergence tests were performed to determine the level of theory and basis set required for accurate 
calculations of three-body hydrogen interactions. We found that the coupled-clusters (CCSD(T)[ p8t) 
level of theory with Dunning's augmented, correlation consistent, triple zeta basis set, aug-cc-pVTZ |29{ 
was sufficient to converge the three-body energies to within 0.1 K. The generalized counterpoise method 
of Valiron and Mayer |Q was used to correct for basis set superposition error (BSSE). 

In order to compare with the three-body contributions calculated using the triple-dipole approximation 
for isotropic oscillators (Eq.4), we have calculated isotropic potentials for both the linear and equilateral 
geometries. The isotropic potential is obtained by weighted averaging of the potentials representing 
each of the ten combinations of orientations. The weight placed on each potential is proportional to 
the degeneracy of the specific orientation. In effect, this corresponds to uniformly averaging over all 27 
possible orientation. This is illustrated for the equilateral geometry in Fig. 4. 

At close range, Hartree-Fock level contributions such as electron exchange and repulsion effects can 
be large ||3l| , |32| , |33| . Only at larger intermolecular separations do dispersion forces dominate. With 
this in mind, we calculated the Hartree-Fock nonadditive three-body potential energy surface using the 
same basis set and BSSE correction. A "dispersion only" three-body potential surface was approximated 
by subtracting the Hartree-Fock PES from the original CCSD(T) surface. Fig. 5 shows the resulting 
dispersion potentials for the equilateral (filled symbols) and linear (empty symbols) geometries. This 
approach is not a true calculation of the dispersion potential. But, the quickly decaying HF potentials 
suggest that the remaining interactions are mainly dispersion. Fitting the long range energies (r'noarost 
> 3.8 A for the linear and r„oarcst > 4.8 A for the equilateral) using the leading triple-dipole dispersion 
function given by Eq. (4) resulted in a root-mean-square error of 0.013. By including the next leading 



term, 

U{DDQ) = ZjDDQ) ^ '°'^^'^ " '^ ^"^^^^^^ + ^ '''f' " ^2) [3 + 5 cos(203)] ^ ^^^^ 

™0 /v.'i /vit: 

'l2'23'31 

one further reduces the error to 0.006. We did not find significant improvements to the fit by including 
the remaining triple dispersion terms. The dispersion model fitted to both the linear and equilateral 
geometries are plotted in Fig. 5 as dashed and solid curves, respectively. The resulting dispersion energy 
constants are Z{DDD) = 49607 K-A^ and Z{DDQ) = 5763 K-A". McDowell and Meath report an 
"exact" value for Z{DDD) of 49833 K-A^, calculated from dipole oscillator strength distribution data p3. 
McDowell [E^ observed that for a several specific orientations of H2 trimers (all three molecules oriented in 
the same plane), nonadditive induction contributions can be as large as those from DDD. This, however, 
is not apparent in the isotropic potentials, where the induction effects have been completely "washed-out" 
by the uniform angle-averaging process [p2[ , as is the case with permanent multipole interactions. 

The Hartree-Fock potentials in Fig. 5 indicate that non-dispersion interactions, such as exchange/repul- 
sion effects, are substantial at short range. We will present a thorough analysis of the different origins of 
hydrogen three-body interactions in a forthcoming article [^ . It is clear from the present study that the 
DDD approximation, alone, is not sufficient at describing the total (CCSD(T)) three-body interactions 
at short distances. Whether this hold true for systems containing the other constituents remains unclear 
and would require separate ah initio calculations. 

4 Discussion 

In this paper we have found that three-body interactions play an important role in altering the interactions 
between adatoms within IC's between carbon nanotubes. The reduction in the well depth of the pair 
interaction is 54%, 28% and 25% for II2, He and Ne, respectively. The magnitude of these shifts is roughly 
similar to the result of a crude estimate, i.e. a factor of three increase over the effect of a single planar 
graphite surface. The DDD term even becomes larger in magnitude than the two-body potential for 
interparticle separations above ~ 7.4 A and 6.3 A for He and H2, respectively. Such a large DDD effect 
implies a significant change in the condensation properties of H2 and He in carbon nanotube bundles. For 
(H2)2 the reduction in well depth is sufficient as to destroy the bound state. The perimeter of a nanotube 
bundle is accessible for adsorbate particles of bigger size. In our approach, the observed reduction in the 
well depth of the pair interaction due to DDD effect is still significant: 35%, 24%, 15%, 19%, 22% and 
28% for H2, He, Ne, Ar, Kr and Xe, respectively. 

The important consequence from the DDD effect for quasi-ld IC phase of H2 molecules and He atoms 
is that there is no dimer bound state. The lack of the latter implies that there is no bound state for the 
N particle system. Then, the ground state of H2 molecules interacting with the screened pair potential 
is a gas, not a liquid. This result raises an interesting question: what is the low temperature behavior of 
such a marginally unbound system? 

Finally, ah initio calculations of three-body interactions for linear and equilateral triangle arrange- 
ments of three H2 molecules provide implications for the behavior of three-dimensional phases of hydrogen. 
Moreover, they justify a future ah initio study of the exchange interaction of two H2 molecules and a 
carbon atom within a nanotube. 

We are grateful to Keith Williams, Paul Sokol, Peter Eklund and Moses Chan for helpful discussions. 
This research has been supported by the Army Research Office under grant DAAD 19 — 99 — 0167. 
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Tables 



Table 1: DDD effect for H2, He and Ne adsorbates confined within IC's of a nanotube bundle. All lengths 
are in A units; e and eg are in meV. 



Quantity 


H2 


He 


Ne 


e{eo) 


1.35(2.96) 


0.67(0.93) 


2.72(3.64) 


5{%) 


54 


28 


25 


Zm(Zam) 


3.55(3.41) 


3.00(2.97) 


3.12(3.09) 


P 


6.3 


7.4 


- 



Table 2: DDD effect for H2, He, Ne, Ar, Kr and Xe adsorbates confined into the groove of a nanotube 
bundle. All lengths are in A units; e and cq are in nieV. The free space potential parameters for Ar, Kr 
and Xe were taken from Ref. 1231. 



Quantity 


H2 


He 


Ne 


Ar 


Kr 


Xe 


a 


2.97 


2.74 


3.00 


3.40 


3.42 


3.40 


e{eo) 


1.93(2.96) 


0.71(0.93) 


3.08(3.64) 


10.1(12.3) 


13.5(17.2) 


17.5(24.2) 


S{%) 


35 


24 


15 


19 


22 


28 


^m\Zom) 


3.49(3.41) 


3.00(2.97) 


3.11(3.09) 


3.79(3.76) 


4.05(4.01) 


4.42(4.36) 



Figure Captions 



(3) 



Fig.l. Free-space potential [/^^ (circle/full curve), the DDD potential V^^^ (triangle/full curve) 

(2'] 

and the effective pair potential F^g (solid curve) for H2 molecules adsorbed within the IC's of nanotube 

bundles are shown. The inset shows the behavior of V^^j close to the sign-reversal point. 

Fig. 2. Same as Fig.l for He. 

Fig. 3. Same as Fig.l for Ne. The inset shows that there is no reversal in sign for V^^^ in this case. 

Fig. 4. Angle-averaging the 27 orientations of the equilateral triangle arrangement of (112)3 to obtain 
an isotropic potential, show in the bold curve with filled circles. The thin curves represent each of the 10 
distinct orientations at each separation. 

Fig. 5 .Hartree-Fock plus triple-dispersion fit of equilateral and linear isotropic potentials. Calcu- 
lated data are represented in symbols (filled = equilateral, empty = linear). Solid curves denote fits to 
equilateral geometry; dashed curves are fits to linear geometry. 
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